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ABSTRACT 

A  focused  system  having  non-coincident  foci  was  used  to  study  the 
phenomenon  of  non-linear  mixing  in  water.   Primary  frequencies  of  343 
kHz  and  383  kHz  provided  a  scattered  difference  frequency  of  40  kHz. 
Enhancement  of  the  difference  frequency  signal  at  incipient  cavitation 
was  verified.   A  delay  between  energization  of  the  system  and  signal 
enhancement  was  linked  to  the  state  of  non-linearity  of  the  medium  and 
provided  strong  evidence  for  the  existance  of  a  virtual  source  at  40 
kHz.   Directivity  measurements  of  the  focused  system  at  40  kHz  revealed 
a  multilobed  structure,  but  the  presence  of  pseudosound  degraded  the 
measurements  and  attempts  to  correlate  the  lobe  structure  with  a  parti- 
cular virtual  source  configuration  were  unsuccessful. 
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I.   INTRODUCTION 

A.    BACKGROUND 

Two  or  more  sound  waves  passing  through  a  common  region  in  a  fluid 
will  experience  an  interaction.   This  coupling  effect  is  due  to  the  non- 
linear terms  in  the  wave  equation  and  one  of  the  results  is  a  scattered 
wave.   The  frequency  of  the  scattered  wave  is  determined  by  the  fre- 
quencies of  the  interacting  waves. 

By  utilizing  Lighthill's  equation  of  aerodynamically  generated 
sound  [1],  Westervelt  was  able  to  mathematically  analyze  the  effect  of 
two  sound  waves  passing  through  a  common  region  in  a  perfect  fluid  [2]. 

When  considering  only  the  lowest  order  scattering  process,  he 
obtained  a  source  function  which  was  quadratic  in  the  primary  field 
variables.   The  solution  of  the  resulting  wave  equation  for  the  plane 
wave  case  showed  that  no  combination  waves  should  be  scattered  out- 
side the  region  transversed  by  two  perfectly  collimated  plane  waves 
forming  a  non-zero  angle  between  them.   Westervelt 's  work  was  corrobor- 
ated when  experiments  by  Bellin  and  Beyer  [3], and  Muir  and  Blue  [4] 
showed  little  deviation  with  his  plane  wave  scattering  theory. 

The  solutions  for  the  interaction  of  two  concentric  cylindrical 
waves  and  of  two  concentric  spherical  waves  was  presented  by  Dean  [5], 
while  the  case  of  eccentric  cylindrical  and  spherical  sources  was 
treated  by  Lauvstad  [6], 

An  experimental  investigation  by  Dunn,  Kuljis  and  Welsby  [7]  was 
directed,  in  part,  towards  studying  the  scattered  sound  from  a  region 
having  pressures  bounding  the  cavitation  region.   To  obtain  high 
acoustic  intensities  with  only  moderate  power,  they  used  a  roughly 


focused  system  comprised  of  two  spherically  concave  transducers  having 
coincident  foci. 

One  transducer  was  driven  at  326  kHz  at  constant  voltage,  while 
the  other  transducer  was  driven  at  373  kHz  and  its  voltage  varied. 
They  found  that  when  the  level  of  fixed  excitation  was  sufficient  to 
cause  cavitation,  the  signal  created  at  the  difference  frequency  (47 
kHz)  rose  smoothly  to  a  plateau  as  the  variable  excitation  was  increased; 
however,  when  the  fixed  excitation  was  insufficient  to  induce  cavita- 
tion, increasing  the  variable  excitation  resulted  in  a  very  slight  in- 
crease in  signal  level  up  to  the  cavitation  threshold.   Beyond  the 
threshold  the  signal  level  rose  rapidly,  attaining  the  same  plateau 
as  in  the  first  case.   One  of  their  conclusions  was  that  the  inter- 
action of  standing  wave  fields  at  two  different  primary  frequencies 
can  produce  a  virtual  source  of  radiated  energy  at  the  difference 
frequency,  and  that  although  cavitation  increases  the  effect,  it  is 
not  necessary. 

B.    OBJECTIVE 

Since  only  a  limited  amount  of  experimental  data  on  focused  systems 
was  available,  a  logical  first  step  in  entering  this  experimental  area 
was  to  confirm  some  of  the  previously  obtained  results.   Once  this  was 
accomplished,  the  goal  was  to  study  some  of  the  phenomena  which  were 
observed  by  earlier  researchers,  but  which  had  received  little  at- 
tention. 

The  work  of  Dunn,  Kuljis  and  Welsby  alluded  to  the  directional 
characteristics  of  the  focused  field,  but  this  aspect  of  a  focused 
system  was  not  studied  in  detail;  therefore,  it  was  decided  to  investi- 
gate the  directional  field  of  a  focused  system  having  non-coincident 
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foci  and  determine  to  what  extent  a  virtual  source  of  the  difference 
frequency  is  created.   Additionally,  some  quantitative  measures  of 
the  time  dependence  of  the  scattered  difference  frequency  was  desired, 
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II.   THEORY  FOR  NON -LINEAR  MIXING 

When  the  wave  equation  is  normally  derived,  approximations  are  made 
which  are  tenable  because  of  the  low  amplitude  signals  usually  en- 
countered in  acoustics  work.   One  example  of  such  an  approximation 
involves  the  Acoustic  Continuity  Equation: 

•— \  s~\  s~\  is  the  instantaneous  total  fluid 

\-s  \0   '     r-^j       density  of  a  small  region 

\0      ~        ambient  density 

V-^|   ~   acoustic  density 

\/      ZT   the  instantaneous  vector  particle 
velocity 

When  equation  (1)  has  been  expanded  the  following  approximations 

may  be  made:  O       _.         /r\ 

-fci  =    -p-  <  <  1 

Po  P* 

Vp.  =    Vp0   -    ^°  =  O 

is  the  instantaneous  total  field 


-*^    y~~.        i^      i—^  J.  o      Lite       luaLaukau^viio      i-wi-ui      i.x^  j.u 

A^  /°     I    '        pressure  at  a  point  in  the  medium 


£)    —   ambient  pressure 


,  zi   acoustic  pressure 


With  these  approximations,  equation  (1)  simplifies  to: 


e-v-v+|B  =  o 


(2) 
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When  this  result  is  coupled  with  the  assumption  of  adiabatic  reversi- 
bility for  an  acoustic  wave,  the  resulting  wave  equation,  equation  (3), 
represents  a  linear  system.   By  linear  is  meant  that  the  system  para- 
meters are  independent  of  the  amplitude  of  excitation;  therefore,  the 
Superposition  Theorem  applies  when  two  or  more  sinusoidal  waves  of  differ- 
ent frequency  exist  in  the  medium. 

T\q    —     bulk  modulus  for  a  static  fluid 
\/   /w  =  velocity   °f  propagation 

When  waves  of  high  intensity  are  generated,  the  previous  approxi- 
mations are  invalid,  and  the  ratio  ^  /  fo  i-n  equation  (3)  should  be 
replaced  by  a  function  '•/P  which  is  a  function  of  acoustic  pressure 
(8). 

*fe    -  ^o0(i+a<p  +  a2p*+...)   <4) 


where  a^,  an,  ....  are  constants 
determined  by  the  parameters  of 
the  fluid 


When  this  substitution  is  made,  the  time-space  relationship  of 
the  pressure  waveform  becomes  dependent  upon  a  velocity  of  propagation 
which  is,  at  any  point,  a  function  of  the  instantaneous  pressure  at  that 
point.   This  non-linear  effect  will  cause  two  sinusoidal  waves  of  fre- 
quencies f,  and  f   to  interact,  generating  frequencies  nf  +  mf 
Here  n  takes  on  all  integral  values  from  0  to  cxO    t    and  m  takes  on  all 
integral  values  from  -  CO  to  +  OO    .   While  this  approach  is  very  general, 

it  expresses  the  concept  of  non-linear  mixing. 
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Westervelt's  precise  approach  was  based  upon  a  calculation  of  the 
scattered  pressure  field  due  to  non-linearities  within  a  small  inter- 
action volume  [2].   Discounting  any  distortion  of  these  second-order 
effects  as  they  propagate  in  a  non-linear  manner,  he  obtained  an  in- 
homogeneous  wave  equation  for  the  acoustic  density  of  the  scattered  wave: 


2 


2-^-'  „Z 


n2es=c;2{D1E-v2^T+ex;2(^  v]}  <» 

I  =  Kinetic  energy  density  =   -r-  T0  LA 
V  =  Potential  energy  density  ■  —Cc   r©  \ 
r~  =  Total  energy  density  ■      -+-  \/ 

Qq   =  Velocity  of  sound 

vv  -   Acoustic  density  of  the  scattered  wave 

[A.  =   Instantaneous  particle  velocity 

The  theoretical  extension  of  Westervelt's  work  to  cylindrical  and 
spherical  spreading  waves  was  accomplished  by  Lauvstad  [6]  and  Dean  [5], 
but  the  theoretical  analysis  of  interactions  within  a  focused,  standing 
wave  system  has,  as  yet,  not  been  accomplished.   Because  of  the  lack  of 
a  theoretical  model,  measurements  and  observations  associated  with  this 
research  could  only  be  compared  to  experimental  work  done  by  others. 
It  was  expected  though,  that  sufficient  observations  would  spark  some 
theoretical  opinions,  expressed  as  extensions  of  the  approaches  dis- 
cussed above. 
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III.   EXPERIMENT  DESCRIPTION 

A.    TRANSDUCERS 

The  transducers  chosen  for  the  focused  system  were  three  inch 
diameter  lead  zirconate-titanate  hemispheres,  having  a  resonant  fre- 
quency of  355  kHz  when  operated  in  the  thickness  made.   The  hemisphere 
was  chosen  in  lieu  of  a  portion  of  a  hemisphere  because  it  was  antici- 
pated that  the  driving  voltage  to  the  system  would  be  limited  by  avail- 
able equipment.   The  additional  surface  area,  therefore,  was  needed  to 
reduce  the  voltage  required  to  achieve  the  desired  cavitation  condition. 
This  decision  immediately  precluded  the  possibility  of  studying  a  focused 
system  with  coincident  foci  since  the  focus  for  a  hemisphere  does  not 
occur  at  an  on-axis  distance  from  the  surface  which  is  greater  than 
the  radius  of  curvature  [9]. 

Mountings  for  the  transducers  were  constructed  from  aluminum  cylin- 
ders, and  the  transducers  were  held  in  place  with  epoxy,  which  also 
acted  as  a  seal.   The  construction  details  are  shown  in  figure  1. 
After  the  transducer  housing  was  completely  assembled,  conducting  silver 
paint  was  applied  around  the  inner  perimeter  of  the  transducer  so  that 
a  continuous  conducting  strip  bridged  the  plastic  insulator  and  contacted 
the  aluminum  housing.   While  the  paint  was  still  tacky,  narrow  strips  of 
aluminum  foil  were  pressed  into  the  paint,  over-lapping  the  edges  of 
the  plastic  insulator.   This  was  done  to  insure  that  an  adequate  cur- 
rent path  from  the  transducer  face  to  the  aluminum  housing  would  be 
formed.   The  aluminum  housing  provided  the  return  or  ground  side  of 
the  circuit.   Finally,  the  complete  unit  was  coated  with  epoxy  and  then 
liquid  neoprene.   Although  the  transducers  had  similar  characteristics 


15 


they  were  numbered  1  and  2  so  that  control  could  be  maintained.   Fur- 
ther reference  to  transducer  1  or  transducer  2  will  be  abbreviated  #1 
or  #2,  and  the  term  "transducer"  will  be  used  to  mean  the  complete  unit 
shown  in  figure  1.   All  distances  from  a  transducer  are  measured  from 
the  mounting  face,  unless  otherwise  specified. 

Two  different  probes  were  used,  one  for  measurements  at  the  rela- 
tively high  primary  frequencies  (343  kHz  and  383  kHz),  and  one  for  the 
relatively  low  frequency  (40  kHz)  difference  wave.   The  active  element 
of  the  high  frequency  probe  was  a  1/8"  square  piece  of  barium  titanate 
mounted  on  a  1/8"  diameter  tube.   The  low  frequency  probe  utilized  an 
Atlantic  Research  LC-65  pressure  transducer  mounted  in  a  brass  housing 
(see  figure  2).   For  the  frequencies  measured,  the  probe  used  was 
smaller  than  a  wavelength  in  water  (  ,A  =  0.392  cm  at  383  kHz  and  ?\   = 
3.77  cm  at  40  kHz).   Both  probes  were  calibrated  by  the  method  of  sub- 
stitution [10]  using  calibrated  LC-32  and  USN  standard  type  E-8  hydro- 
phones.  The  sensitivities  of  the  transducers  were  as  follows: 

DIFFERENCE,  OR 

LOW  FREQUENCY  PROBE      -91  dB  re  lV//ibar  at  40  kHz 

PRIMARY,  OR  -102  dB  re  lV/yubar  at  343  kHz 

HIGH  FREQUENCY  PROBE     -96.2  dB  re  lV/yubar  at  383  kHz 

The  directivity  plot  for  the  low  frequency  probe,  figure  3,  shows 
that  the  beamwidth,  based  upon  the  -3  dB  points,  is  at  least  30°. 
The  high  frequency  probe  was  essentially  omnidirectional, 

B.    ELECTRONICS 

The  following  types  of  equipment  were  used  in  the  investigation: 

Item  Model 

SIGNAL  GENERATOR  HP  650A 

FREQUENCY  COUNTER  HP  5232A 


16 


OSCILLOSCOPE  TEKTRONIX  515A 

VOLTAGE  AMPLIFIER  HP  450A 

POWER  AMPLIFIER  GR  1233A 

WAVE  ANALYZER  GR  1900A 

*WAVE  ANALYZER  HP  3590A 

RADIO  TRANSMITTER  TDE-3  (navy  model) 

VTVM  HP  400C 

X  -  Y  PLOTTER  HONEYWELL  320 

*was  not  available  until  late  in  the  investigation 

The  TDE  -3  is  a  military  transmitter  capable  of  producing  125  watts 
in  the  CW  mode.   Because  it  was  designed  for  communications  purposes, 
it  was  not  possible  to  drive  a  transducer  for  longer  than  30  seconds 
without  risking  damage  to  the  transmitter.   This  restriction  limited 
the  use  of  the  transmitter  to  investigations  requiring  short  observa- 
tion. 

Measurements  made  with  the  GR  1900A  were  over  a  50  Hz  bandwidth, 

and  the  meter  response  time  was  set  at  either  fast  (0.15  second),  or 

# 
medium  (0.5  second).   The  997<,  confidence  limits   for  the  meter  speeds 

are : 

Fast  -2.5  to  +3.3  dB 

Medium  -1.4  to  +1.7  dB 

Because  the  signal  generators  used  had  a  tendency  to  drift,  all 
measurements  were  made  using  AFC. 


One  chance  in  100  that  the  ratio  of  the  long  time  average  value 
to  a  selected  meter  reading  will  be  beyond  the  range  shown. 
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C.    MECHANICAL 

A  set  of  rails  with  two  moveable  slides  was  used  to  position  the 
probes.  A  special  holder,  constructed  in  two  halves,  allowed  the  probe* 
or  transducers  to  be  positioned  and  then  rotated  to  fixed  angular  posi- 
tions.  Figure  5  shows  one  transducer  and  the  low  frequency  probe  in^ 
stalled  on  the  rail. 

To  investigate  the  directional  properties  of  the  focused  system, 
it  was  considered  easier  to  rotate  the  focused  system  than  to  move  the 
probe  to  various  angular  positions.  A  jig  was  constructed  which 
allowed  the  focused  system  to  be  rotated  approximately  90°  for  a  7  cm 
transducer  separation.   Figure  6  shows  the  focused  system  in  the  jig. 
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IV.   EXPERIMENTAL  PROCEDURE;  RESULTS 

A.  GENERAL 

Originally  the  investigation  was  to  have  been  carried  out  in  a 
2'  x  4'  x  6'  plywood  tank  lined  with  one  inch  of  rubberized  hair;  but, 
because  initial  measurements  showed  that  standing  waves  were  present, 
a  61  x  8'  x  24 '  anechoic  tank  was  used  instead. 

The  mechanical  resonance  of  each  transducer  was  determined  by  find- 
ing the  frequency  which  maximized  the  output  under  the  condition  of  con- 
stant voltage  input  [11].   This  value  was  found  to  be  355  kHz  for  each 
transducer.   Equipment  limitations  prevented  analyzing  a  difference 
frequency  of  greater  than  50  kHz,  and  so  a  40  kHz  difference  was  de- 
cided upon.   The  choice  of  the  primary  frequencies  was  made  by  examining 
the  loaded  output  of  the  TDE  -  3  for  various  frequencies  and  noting  those 
which  minimized  the  VSWR  on  the  connecting  coaxial  cable.   The  most  ac- 
ceptable VSWR,  1.1:  1,  was  realized  at  343  kHz  and  so  the  TDE  -  3  was 
used  at  this  frequency  while  the  drive  at  383  kHz  was  achieved  through 
the  arrangement  shown  in  figure  4. 

Because  this  investigation  was  concerned  with  the  interaction  of 
acoustic  fields,  it  was  important  to  know  the  directional  characteris- 
tics of  each  transducer. 

B.  TRANSDUCER  DIRECTIVITY 

Transducer  #1  was  driven  with  a  15V  p-p,  383  kHz  signal,  and  the 
high  frequency  probe  was  located  at  5.4  cm  from  the  transducer.   The 
transducer  was  then  rotated  about  its  support,  and  the  directional  pat- 
tern shown  in  figure  7  was  obtained.   A  second  directional  pattern  was 
obtained  at  343  kHz  with  the  probe  at  10.8  cm  from  the  transducer  (see 
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figure  8).  An  attempt  to  repeat  the  directional  measurements  at  drive 
levels  producing  cavitation  was  unsuccessful  because  of  large  noise  - 
like  fluctuations. 

C.  ON-AXIS  INTENSITY 

O'Neil's  theory  predicting  the  on-axis  intensity  for  a  spherical 
focusing  radiator  is  not  valid  for  a  hemisphere  because  one  requisite 
is  that  the  radius  of  the  circular  boundary  of  the  spherical  cap  be 
large  relative  to  the  depth  of  the  concave  surface  [9];  therefore,  it 
was  decided  to  probe  the  on-axis  near-field  to  get  a  rough  estimate  of 
the  relative  intensity.   The  on-axis  plots  were  similar  for  both  pri- 
mary frequencies,  and  figure  9  shows  the  plot  for  #2  at  383  kHz  as  a 
representative  result.   For  the  same  driving  voltage,  the  level  of  the 
signal  at  343  kHz  was  greater  than  at  383  kHz.   This  result  was  expected 
though,  since  operation  near  resonance  is  more  efficient. 

D.  LOCATION  OF  FOCUS 

Theory  predicts  that,  for  a  hemisphere,  the  focus  will  be  located  at 
the  same  point  as  the  center  of  curvature.  This  was  verified  by  placing 
transducer  #1  in  a  14  cm  x  19  cm  x  28  cm  plexiglass  tank  and  then  driving 
it  at  343  kHz  with  a  170V  p-p  signal.   This  level  of  drive  was  sufficient 
to  achieve  transient  cavitation  [12].  Attendant  with  this  effect  was 
sporadic  acoustic  streaming  which  carried  cavitation  bubbles  through  the 
focal  region,  allowing  visual  observation  of  both  the  convergence  and 
the  divergence  around  the  focus.   The  effect  was  intermittent,  and  hence 
it  was  difficult  to  record  on  a  photograph.   Figure  10  shows  the  micro- 
bubbles  in  the  focal  region  just  before  foggy  streaming  [13]  commenced. 
A  measurement  of  the  on-axis  location  of  the  focus  confirmed  the  fact 
that  it  was  within  0.1  cm  of  the  center  of  curvature  of  the  hemisphere. 
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E.    CAVITATION  THRESHOLD 

1.    Visual  Determination 

Because  of  the  presence  of  small  cavitation  bubbles  between 
the  focused  transducers  enhances  the  level  of  the  scattered  difference 
frequency,  it  was  necessary  to  have  a  rough  estimate  of  the  driving 
voltage  required  to  achieve  cavitation.   This  was  accomplished  by  focusing 
the  two  transducers  and  then  placing  the  focused  system  in  the  plastic 
water  tank.   Future  references  to  the  focused  system  shall  be  interpreted 
to  mean  the  two  transducers  with  parallel  faces,  separated  by  7.0  cm. 
For  the  test,  #1  was  operated  at  343  kHz  and  #2  at  383  kHz.   The  driving 
voltage  applied  to  #2  was  kept  at  zero  while  the  driving  voltage  to  #1 
was  increased  from  zero  in  5V  steps.  At  35V  p-p,  small  bubbles  were  seen 
in  the  vicinity  of  the  focus  of  #1  and  after  a  5  minute  wait  bubbles  were 
found  to  be  fairly  uniformly  distributed  between  the  two  transducers. 
Some  of  the  bubbles  were  seen  to  vibrate,  both  horizontally  and  vertically. 
While  the  predominant  motion  was  vertical  (due  to  buoyant  forces),  some 
short  horizontal  movements  were  witnessed,  both  with  and  without  coalesc- 
ence.  Of  particular  interest  was  the  observation  that  within  the  space 
separating  the  two  transducers,  the  vertical  bubble  motion  occurred  in 
discrete  jumps  between  visually  undefined  layers. 

When  the  voltage  to  #1  was  reduced  to  zero,  it  was  found  that 
incipient  cavitation  occurred  when  the  voltage  to  #2  was  30V  p-p.   Fig- 
ure 11  shows  some  of  the  bubbles  in  the  vicinity  of  the  focus  of  #2  one 
minute  after  driving  voltage  was  applied. 

To  ascertain  whether  the  voltage  applied  to  one  of  the  transducers 
in  the  focused  system  would  influence  the  cavitation-voltage  threshold 
of  the  other  transducer,  various  combinations  of  driving  voltage  were 
tried.   The  following  table  summarizes  the  results  and  shows  that  the 
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cavitation-voltage  threshold  for  one  transducer  in  a  focused  system  is 
strongly  influenced  by  the  voltage  applied  to  the  other  transducer. 


Transducer  #1 

Trat 

lsducer  #2 

Cavitation 

Cavitation 

@  343  kHz 

@ 

383 

kHz 

at  focus 

at  focus 

Volts  p-p 

Volts 

P-P 

of  #1 
no 

of  #2 

30 

0 

no 

35 

0 

yes 

no 

0 

25 

no 

no 

0 

35 

no 

yes 

10 

30 

yes 

yes 

15 

30 

yes 

yes 

25 

30 

yes 

yes 

25 

4 

no 

no 

25 

15 

no 

yes 

2.    Subharmonic  Detection 

Late  in  the  conduction  of  the  investigation  a  wave  analyzer  with 
the  capability  to  analyze  high  frequencies  was  received  and  so  it  was 
used  to  recheck  the  cavitation-voltage  threshold  by  the  method  of  sub- 
harmonic  detection  [14].   The  system  was  focused  and  the  high  frequency 
probe  was  placed  at  the  focus  of  #1.   The  probe  output  went  directly  to 
the  wave  analyzer  which  was  tuned  to  one-half  of  the  driving  frequency 
(1st  subharmonic)  of  #1,   No  voltage  was  applied  to  #2  during  this  part 
of  the  test.   The  drive  to  #1  was  increased  until  the  subharmonic  was 
detected.   This  procedure  was  carried  out  for  #1  and  #2  at  both  primary 
frequencies  with  the  following  results: 

Cavitation-voltage    Cavitation-voltage 
threshold  at  threshold  at 

343  kHz  383  kHz 

Transducer  1  18V  p-p  22V  p-p 
Transducer  2                 14V  p-p  18V  p-p 
Subharmonic  generation  is  acknowledged  to  be  the  most  accurate  indi- 
cation of  incipient  cavitation  because  bubbles,  if  visible  at  all,  do  not 
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occur  immediately  when  the  threshold  is  reached.   For  this  reason  the  values 
obtained  by  subharmonic  detection  were  accepted  as  the  voltage  threshold. 

F.    SPL  AT  FOCUS 

An  attempt  was  made  to  measure  the  sound  pressure  level  (SPL)  of 
the  primary  frequencies  at  the  focus  of  each  transducer  while  the  system 
was  focused.   Figure  12  shows  the  plot  of  SPL  vs  Driving  Voltage  for  each 
transducer  at  each  primary  frequency.   The  accuracy  of  these  values  is 
questionable  since  the  presence  of  the  probe  alters  the  field  configura- 
tion and  eliminates  the  central  portion  of  the  focus;  nevertheless,  the 
curves  provide  some  useful  information  since  the  shapes  still  reveal  the 
relationship  between  driving  voltage  and  SPL. 

G.    DIFFERENCE  FREQUENCY  ENHANCEMENT 

Prior  to  the  actual  investigation  of  the  directional  pattern  of  the 
focused  system,  it  was  decided  to  measure  the  SPL  of  the  difference  fre- 
quency to  confirm  that  the  system  actually  did  behave  as  reported  by  Dunn, 
Kuljis  and  Welsby  [7].   The  transducers  were  focused  and  the  probe  was 
placed  at  11  cm  along  the  perpendicular  bisector  of  the  system  axis 
(defined  by  a  line  connecting  the  centers  of  the  two  hemispheres).   The 
driving  voltage  to  #1  was  varied  from  4V  p-p  to  50V  p-p  while  the  voltage 
to  #2  was  kept  constant.   Five  such  runs  were  made  using  constant  values 
of  4,  8,  15,  25  and  30V  p-p.   The  results,  shown  in  figure  13,  indicate 
that  signal  enhancement  does  occur  when  the  cavitation  region  is  reached. 
Figure  14  shows  the  form  of  the  received  signal  prior  to  filtering. 

H.    SYSTEM  DIRECTIONAL  CHARACTERISTICS 

The  directional  characteristics  of  the  scattered  difference  frequency 
were  investigated  by  mounting  the  focused  system  in  the  rotatable  jig 
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(figure  6)  and  making  measurements  at  a  distance  of  66  cm.   The  total 
angular  rotation  of  the  jig  was  limited  to  about  90°  because  of  inter- 
ference between  the  pivot  support  bar  and  the  transducer  support  stems. 
Figures  15  -  18  show  the  radiation  pattern  over  the  range  318°  to  48°  for 
various  combinations  of  driving  voltage.   The  appearance  of  incomplete 
large  lobes  at  the  extremes  of  the  angular  rotation  of  the  jig  led  to  an 
attempt  to  complete  the  lobe  data.   Unfortunately,  when  the  system  was 
readjusted  to  study  the  lobe  at  048°,  the  orientation  of  the  individual 
transducers  shifted  and  the  original  pattern  could  not  be  reconstructed 
(see  figure  18).   One  puzzling  outgrowth  of  the  attempt  was  the  occur- 
rence of  a  portion  of  a  major  lobe  at  090°. 

I.    PSEUDOSOUND  EFFECTS 

The  possibility  that  pseudosound  (non-linear  mixing  of  the  primary 
frequencies  at  the  probe)  was  being  mistaken  for  the  scattered  differ- 
ence frequency  was  considered.   One  brief  experiment  with  SOAB  (a  B.  F. 
Goodrich  Co.  rubber  having  selective  frequency  absorptive  properties) 
to  determine  the  source  of  the  difference  frequency  was  unsuccessful  be- 
cause the  attenuation  of  the  primary  and  difference  frequencies  were 
similar.   A  second  attempt  to  isolate  the  source  of  the  40  kHz  signal  in- 
volved comparing  the  strength  of  the  difference  frequency  while  the 
system  was  focused  with  that  measured  when  both  transducers  were  pointing 
at  the  probe.   If  pseudosound  was  the  source  of  the  difference  frequency, 
then  the  latter  situation  should  have  maximized  the  40  kHz  signal.   The 
results  tended  to  discount  the  effects  of  psuedosound  since  the  received 
signal  for  the  focused  condition  produced  signals  which  were  from  2  to 
20  dB  greater  than  for  the  unfocused  configuration  (the  larger  driving 
voltages  produced  larger  differences).   Trying  another  approach,  the  HP 
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359A  wave  analyzer  was  used  to  obtain  the  focused  system's  directivity 
pattern  at  each  primary  frequency  while  the  system  was  being  driven 
(figures  20  and  21).   The  peaks  and  dips  on  the  plot  were  interpreted  as 
an  interference  effect  due  to  direct  radiated  waves  from  a  transducer 
combining  (at  the  probe)  with  waves  reflected  from  the  other  transducer. 
The  results  only  tended  to  strengthen  the  case  for  psudosound  since  the 
SPL's  of  the  primary  frequencies  were  generally  above  the  SPL  of  the  differ- 
ence frequency.   For  example,  calculations  for  the  conditions  existing 
when  figure  19  was  plotted  show  that  the  following  readings  occur  within 
+2°  of  090°. 

SPL  =  10 . 2  dB  for  383  kHz  at  90° 
SPL  =  30.5  dB  for  343  kHz  at  88° 
SPL  -  0,0  dB  for  40  kHz  at  92° 

Having  found  no  way  to  discern  between  scattered  sound  and  psuedo- 
sound,  attention  was  returned  to  interpreting  figures  15  through  18. 
The  large  lobes  at  318°  and  048°  strongly  suggested  that  the  wave  inter- 
action was  taking  place  primarily  at  the  foci,  creating  two  virtual 
sources.   The  system  geometry  seemed  to  suggest  that  partial  shielding 
of  these  sources  was  influencing  the  directivity.   Figure  22  shows  the 
system  geometry  and  includes  some  limiting  ray  paths,  assuming  point 
sources  at  the  two  foci.   While  it  appeared, at  first,  that  the  measured 
directivity  patterns  could  be  explained  by  two  virtual  sources  at  the 
foci,  the  explanation  was  clouded  by  many  unknowns.   The  strength  and 
phase  relationship  of  the  virtual  sources  was  not  known.   More  important, 
the  physical  size  of  the  virtual  sources  was  unknown. 

J.    TIME  DEPENDENT  EFFECTS 

The  time  dependence  of  SPL  measurements  was  obvious  throughout  the 
experiment.   Although  it  was  easy  to  rationalize  the  variation  of 

25 


measurements  over  long  periods  of  time  as  being  caused  by  reflections  from 
large  bubbles,  the  short  term  fluctuations  required  further  investigation. 

The  low  frequency  probe  was  placed  14  cm  from  the  center  of  the  focused 
system,  at  an  angular  position  of  318°  and  was  aimed  at  the  focus  of  trans- 
ducer #1.   The  probe  output  was  filtered,  and  the  40  kHz  component  was 
converted  to  d.c.  to  drive  the  x-y  plotter,  thus  recording  the  differ- 
ence frequency  voltage  level  vs.  time.   Figures  23  to  25  show  the   SPL 
versus  time,   for  four  combinations  of  driving  voltage.   Figure  23  indi- 
cates that  there  is  a  definite  time  delay  between  the  application  of  the 
system  voltages  and  the  enhancement  of  the  difference  frequency.   Figures 
24  and  25  reveal  that  higher  driving  voltages  eliminate  the  delay  time 
but  seem  to  produce  an  overshoot  of  the  final  average  value.   Visual 
observation  did  not  allow  complete  correlation  of  signal  level  and  bub- 
ble activity,  but  the  absence  of  any  visible  bubbles  during  the  first 
couple  of  seconds  after  turn-on  seem  to  allow  the  following  interpreta- 
tions of  figure  24. 

The  system  is  energized  at  time  t   and  cavitation  nuclei  begin  grow- 
ing, thus  increasing  the  compressibility  of  the  water.   This  non-linear 
action  increases  until  time  t. ,  when  bubble  growth  and  coalescence  begin 
to  outweigh  the  creation  of  new  cavities.   As  this  process  continues,  the 
radiation  resistance  of  the  medium  decreases  and  the  radiated  power  is 
reduced  [15].   At  time  t.  a  fairly  stable  condition  is  reached  with 
signal  fluctuation  being  caused  by: 

1)  Interference  effects,  due  to  energy  scattered  by  the  bub- 
bles combining  out  of  phase  with  the  virtual  sources. 

2)  Bubbles  screening  a  portion  of  the  interaction  region, 
thereby  altering  the  virtual  sources. 

3)  The  absorption  of  energy  by  the  bubbles. 
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V.   CONCLUSIONS 

The  directivity  plot  for  the  difference  frequency  in  a  focused  system 
was  found  to  have  a  multilobed  structure.   While  the  magnitude  of  the 
driving  voltages  altered  the  angular  position  of  the  lobes  very  little, 
the  size  of  the  lobes  changed  unpredictably.   Increasing  the  distance  at 
which  measurements  were  made  shifted  the  lobes  a  few  degrees,  indicating 
that  interference  effects  might  be  responsible  for  the  lobe  structure. 

Attempts  to  relate  the  lobe  structure  to  some  virtual  source  con- 
figuration in  the  interaction  region  were  inconclusive.  Although  it  was 
suspected  that  the  virtual  source  would  occupy  a  large  portion  of  the 
volume  between  the  transducers,  it  seemed  likely  that  the  energy  density 
of  the  source  would  be  greatest  around  the  two  foci.   If  this  were  true, 
then  the  shape  of  the  radiation  pattern  should  have  been  dominated  by  the 
effect  of  these  "point1"  sources.   Furthermore,  ray  theory  showed  that  the 
major  lobes  should  appear  somewhere  between  322°  and  038°.   Probing  of 
the  field,  however,  showed  the  major  lobe  in  the  first  quadrant  occurred 
at  about  090°,  a  location  where  the  difference  frequency  signal  was  ex- 
pected to  be  smallest.   A  follow-up  comparison  of  the  difference  frequency 
SPL  with  that  of  the  primary  frequencies  over  the  positions  from  15°  to 
95°  indicated  that  the  SPL  at  the  difference  frequency  was  always  10  to 
35  dB  lower  than  the  SPL  for  either  primary  frequency.   With  such  large 
differences,  pseudosound  was  a  possibility  which  could  not  be  dismissed. 

In  spite  of  the  possible  effects  of  psuedosound,  the  case  for  a 
dominant  "point"  source  was  supported  by  the  evidence  that  for  a  given 
position  the  response  of  the  probe  to  the  difference  frequency  could  be 
maximized  by  aiming  the  probe  face  at  the  focus  of  the  transducer  with  the 
largest  driving  voltage. 
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Future  experiments  should  consider  the  presence  of  pseudosound  and 
utilize  techniques  to  circumvent  or  isolate  its  effects.   One  method  might 
be  to  make  directivity  measurements  using  two  probes  having  different  non- 
linear characteristics  and  then  use  a  correlation  technique  to  separate 
the  pseudosound  from  the  scattered  sound. 

The  relative  positions  of  the  two  transducers  in  the  focused  system 
are  critical  and  therefore  apparatus  which  allows  acoustic  measurements 
to  be  made  without  disturbing  the  system  is  needed.   Specifically,  the 
directional  characteristics  for  360°  around  the  system  should  be  obtain- 
able without  moving  the  focused  system.   This  means  that  accurate  posi- 
tioning of  a  probe  around  the  system  is  necessary.  Additionally,  the 
relative  positions  of  the  transducers  should  be  controllable  to  about 
+lmm.   Such  accuracy  would  allow  an  analysis  of  the  effect  of  transducer 
spacing  and  position  upon  the  scattered  difference  frequency. 

The  repeatability  of  acoustic  measurements  made  during  cavitation 
is  very  poor  because  the  cavitation  bubbles  mask,  scatter,  and  absorb 
energy  in  a  random  manner. 

Time  recordings  of  the  difference  frequency  signal  show  that  the 
mechanism  responsible  for  signal  enhancement  is  not  activated  immediately 
after  the  system  is  energized.   Using  the  drive  combination  of  4V  p-p 
and  30V  p-p  at  the  primary  frequencies,  enhancement  was  not  effective 
until  8  seconds  after  turn  on.   The  delay  was  reduced  to  about  1  second 
when  the  driving  combination  was  increased  to  20V  p-p  and  75V  p-p.   This 
phenomenon  invites  a  study  of  the  time  history  of  cavitation  bubbles, 
using  a  high  speed  sequence  camera.   Such  a  study  could  shed  more  light 
on  the  precise  conditions  required  for  signal  enhancement. 
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The  proof,  then,  that  non-linear  mixing  was  actually  achieved  within 
the  medium  is  supported  by  the  fact  that  it  was  possible  to  influence  the 
generation  of  non-linear  conditions  in  the  medium.   By  controlling  the  time- 
rate-  of-  growth  of  the  non-linearity,  the  enhancement  of  the  difference 
frequency  was  similarly  controlled.   Positive  proof  of  non-linear  mixing 
will  exist  only  after  further  research  confirms  that  large  increases  in 
the  SPL  of  the  primary  frequencies  do  not  occur  concurrently  with  the 
creation  of  the  non-linearity. 

While  the  focused  system  currently  has  no  practical  applications,  it 
is  invaluable  as  a  research  tool  for  studying  non-linear  acoustic  phenom- 
ena as  well  as  cavitation  effects.   Perhaps  future  work  can  incorporate 
the  use  of  optical  diffraction  to  determine  the  configuration  of  the  virt- 
ual source. 

Lastly,  a  measure  of  the  efficiency  of  the  virtual  source  is  needed. 
If  it  turns  out  to  be  as  inefficient  as  for  plane  wave  mixing  [16],  then 
the  application  of  a  focused  system  may  depend  solely  upon  its  directional 
properties. 
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APPENDIX  A:   FIGURES 
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Figure  5 

One  Focusing 
Transducer  and  the 
Low  Frequency  Probe 
Mounted  on  the  Rail 


Figure  6 

Focused  System 
Mounted  in  Rotation 
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Figure  10 

Transient  cavitation  in  the  Focal  Region, 
is  170V  p-p  at  343  kHz 

Transducer  at  Right  was  Open  Circuited 
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Figure  14 

Typical  Received  Signal  Prior  to  Filtering 
The  Carrier  is  363  kHz 
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